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We associate to every quantum channel T acting on a Hilbert space Ti a pair of Hamiltonian operators over 
the symmetric subspace of 7i® 2 . The expectation values of these Hamiltonians over symmetric product states 
give either the purity or the pure state fidelity of T. This allows us to analytically compute these measures for a 
wide class of channels, and to identify states that are optimal with respect to these measures. 
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I. INTRODUCTION 

The study of open quantum systems (l| is of interest in 
fields as diverse as quantum information science 0], quan- 
tum control 0|, and foundations of quantum physics |4j]. Let 
T e CP(W) be a completely positive trace-preserving quan- 
tum map i.e., a channel over the finite-dimensional quantum 
state-space fi. The channel T has a (non-unique) Rraus oper- 
ator sum representation 1 5] 

T(X) = J2 A i XA t l^eEnd(ff)] (1) 

i 

where the Rraus operators Ai satisfy the constraint 
^\ A\A - l = 1, which guarantees preservation of the trace 
of a state (density operator) X = p. A fundamental prop- 
erty of a state is its purity p[p] = Tr(p 2 ). States are called 
pure iff p = 1 and mixed if p < 1. In the paradigmatic 
scenario of open quantum systems, a state starts out as pure, 
p = \tp) (ip\, and is then mapped, e.g., via the interaction with 
an environment, to a mixed state by the action of a channel T: 
p[T(p)] = Tt[T 2 (\iP)(tJj\)] < 1. In this case we say that the 
state \ip) has been decohered by the channel. A typical goal 
of, e.g., quantum information processing, is to maximize the 
purity of a state that is transmitted via some channel T. To this 
end a variety of decoherence-reduction techniques have been 
developed, such as quantum error correcting codes (QECCs) 
000E3] an d decoherence-free subspaces (DFSs) Ell- 
in this work we are interested in the intrinsic purity of quan- 
tum channels: 

Definition 1 The purity of the channel T over the subspace 
CcWis 

P(T,C) := min TY[T 2 (|V>XV|)]. ( 2 ) 

Minimization is required since we must consider the worst- 
case scenario. We invoke subspaces in our definition since we 
know from the theory of QECCs and DFSs that it is possible 
to encode quantum information in a manner that maximizes 
purity by restricting to a subspace. In particular: 
Definition 2 If P(T,C) = 1 we say that that C is a 
decoherence-free subspace with respect to T, in short a T- 
DFS. 



In many cases it will not be possible to find a T-DFS. A 
central question we shall be concerned with here is the char- 
acterization of those states that optimally approximate a T- 
DFS, i.e., those states for which P(T, C) is as large as possi- 
ble. Thus: 

Definition 3 The optimal purity of T is 

P(T) := max P(T,C) (3) 

CCH 

Note that P(T) = 1 <^> the set of T-DFSs is non-empty. 
However, this situation is rather rare and generally requires 
that there be a symmetry in the system-environment interac- 
tion. Associated with this symmetry is a conserved quantity: 
quantum coherence. This in turn leads to the preservation of 
quantum information. Here we wish to depart from the no- 
tion of a strict symmetry and explicitly consider the situation 
where one can only expect optimal, as opposed to ideal pu- 
rity. However, the optimization problem defined by P(T) is a 
hard one, since it involves a search over all possible subspaces 
C C H; the number of such subspaces grows combinatorially 
in the dimension of H, which itself may be exponential in the 
number of particles, in a typical quantum information process- 
ing application. Moreover, even if one restricts the problem to 
the computation of P(T, C) (for a given, fixed subspace), one 
is still faced with a complicated-looking functional. 

In this work we focus on the the computation of P(T,C) 
and we associate a Hamiltonian to each channel. This "chan- 
nel Hamiltonian" is a mathematical trick, rather than a physi- 
cal Hamiltonian. But, as we shall show, this has the advantage 
that it allows us to cast the purity problem into the familiar 
framework of computing eigenvalues of Hermitian operators. 
In addition, we show that our channel Hamiltonian leads to 
an elegant physical (re-) interpretation of the channel purity 
in terms of the expectation value of the SWAP operator. 

Our work is also related to questions about channel capac- 
ity; indeed recently it has been shown that multiplicativity of 
generalized maximal purities implies additivity of the mini- 
mal output entropy of the quantum channel. The latter in turn 
is equivalent to the additivity of the Holevo channel capacity 
10. 

We introduce the first channel Hamiltonian in Section ITT1 
We then derive a number of properties and bounds on the pu- 
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rity based on this formalism in Section [TO] We then devote 
Section lTVl to a number of examples designed to illustrate our 
formalism, and derive some interesting properties for a class 
of channels. In Section [V] we derive an alternative interpre- 
tation of the expression for the channel purity, in terms of a 
dual map. It turns out that the same methods we introduce 
for the channel purity also apply to the pure state fidelity of 
the channel. In particular, we can introduce a second channel 
Hamiltonian to this end. This is addressed in Section lVll We 
conclude in Section rvTTl 



II. A HAMILTONIAN OPERATOR FOR QUANTUM 
CHANNELS 



Associated to the channel T we define an operator over 
Ti® 2 : 

Definition 4 The channel purity-Hamiltonian is: 

fi(r) := Yj n h n a ' n y = A A i ■ ( 4 ) 

ij 

(We shall refer to f2(T) simply as the "channel Hamilto- 
nian" until our discussion of the pure state fidelity in Sec- 
tion lVIl ) It follows immediately from fij^. = flji that f2(T) is 
a symmetric, Hermitian operator. Thus f2(T) has the status of 
a Hamiltonian over Ti® 2 . Moreover, f2(T) is independent of 
the particular Kraus operator sum representation chosen for T: 
all possible operator sum-representations of T are obtained by 
considering new Kraus operators of the form A\ = J2j UijAj 
where the Uij's are the entries of unitary matrix. By inserting 
this expression into the definition @ one can explicitly check 
that fi(T) is invariant. 

We now come to our key result: a representation of the pu- 
rity of quantum channels as the expectation value of the chan- 
nel Hamiltonian: 

Proposition For every quantum channel T and subspace C 
one has the identity 

P(T,C) = min(V>® 2 \n(T) I ip® 2 ). (5) 

Wee 

Proof. One has 

TrT 2 (|V^|) = J2(M A M)W A t A M 

ij 

= J2m A } A m 2 = TrU)^ 2 ^ A\A 3 ] 

ij 

= Ti-U)(^ 2 n(T)}. (6) 

Equation (|5} now follows by taking the minimum over S 
Cm 

Note that f2(T) is a formal Hamiltonian over the "double" 
Hilbert space Ti® 2 , and is therefore unrelated to the physical 
Hamiltonian for the original problem. However, as we show 
below, there does exist an attractive physical interpretation of 
Eq.[3] in terms of the expected value of the SWAP operator. 



III. BOUNDS AND OTHER CHANNEL PROPERTIES 

We now derive upper and lower bounds on the purity and 
then give a characterization of T-DFSs. 
Proposition 1 Let uiq (A) denote the minimum eigenvalue of 
the symmetric operator A in the symmetric subspace of Ti. ® 2 , 
and let n + (C) denote the normalized projector over the sym- 
metric part of of C® 2 . Then the following bounds hold: 

Tr[II + (C)O(T)] > P(T,C) > uj+[n(T)}. (7) 

Proof. Note that since O(T) is a symmetric operator, the 
symmetric subspace of Ti® 2 is f2(T)-invariant. Therefore the 
minimum expectation value of 0(T) in this subspace coin- 
cides with the minimum eigenvalue u>q , The lower bound in 
Eq. @ is simply due to the fact that minimization over the 
symmetric subspace of Ti® 2 includes the minimization over 
the product states \ip)® 2 e C® 2 . The upper bound in Eq. Q 
derives from the identity J c \i^)('4>\® 2 = H + (C) (integration 
over the uniform distribution over C 0) and from the obvi- 
ous fact that the average value of a function is no smaller than 
its minimum value. ■ 

Lemma 1 Let T be unital [T(l) = 1]. Then V|V>) G Ti we 
have: 

||fXT)|^ 2 || < 1- (8) 

Proof. Let := ||Atj4j|^i)||. One has the following nor- 
malization condition 

5>l = j2(M A J A ] A m = i>i4ai^>) = i, 

ij ij i 

where in the first (second) equality we used the unitality (CP- 
map) condition J^j A i A \ = 11 (Hi A \ A i = Now 

\MTM® 2 \\ = ||]>>}^)®4a,#}|| 

ij 

< ^114^1^11114^1^)11=5^ 

ij ij 

<$>| = 1, (9) 

ij 

where in the last line we used the Cauchy-Schwartz inequality 
for the Hilbert-Schmidt product of matrices, 

J^PijPji = TrP 2 = (P,pt ) <|| F ||||pt|| 

ij 

ij 

■ 

We now proceed to characterize T-DFSs. To this end we 
introduce a special subspace: 

Definition 5 The subspace Ti n of 17-invariant states (Ti n C 
Ti® 2 ) is the eigenspace of £1 with eigenvalue one. 
Proposition 2 
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i) If V|V>) G C and Vt it holds that: A^ifj) = 
aiU\ip), A\\ip) = a*U^\tp), where U is unitary, then 
|V) 02 G H n . 

ii) Let T be unital. Then C is a T-DFS C® 2 C ft n . 
iii) T-DFS the first inequality in Eq. Q is an equality. 

Proof, i) Notice first that from the CP-map condition, 

J2i A t A i = !. k follows that £V M 2 = l - Now for 

\tp) G C one has that (^ 2 |«(T)|i/;® 2 ) = J2i 3 \ a i a j\ 2 = 

(£,M a ) 2 = i- 

ii) (=>) If C is a T-DFS thenmiri| V;)6C (V'® 2 |fi(T)|^® 2 ) = 
1. But from the Cauchy-Schwartz inequality and Lemma 1 
above one has that 2 \fl(T)\i>® 2 ) < 1 (V|V)) and 
the equality holds iff f)(T)|-0® 2 ) = \^® 2 ) (VjV) G C). 
Now, if |$) is in the symmetric part of C® 2 , one has that 
|*> =II + (C)|tf) = «(C)/ C l^^ 2 )^® 2 !*) [where a(C) := 
dimC(dimC + 1)]; therefore 

fi(T)|¥) = a{C) [ n(T)\Tp® 2 )(?p® 2 \y) 
Jc 

= a(C) J c \^ 2 }(iP® 2 \V). (10) 

It follows that I*) e H n . 

(<J=)IfV|^) eCit holds that 1 1/')® 2 G Hf, then C is clearly 
a T-DFS, i.e., P(T,C) = 1. A fortiori this holds if all the 
elements of the symmetric part of C ® 2 are in TL . 

iii) We have just seen that VI*) = \ip)® 2 such that e C 
one has (^^(T)!*) = 1. By integrating over \ip) one obtains 
that the average [leftmost part of Eq. @] coincides with the 
minimum [middle term in Q]. ■ 



IV. EXAMPLES 

We now present a variety of examples to illustrate our for- 
malism, to actually compute the purity of a number of inter- 
esting channels, and to find the corresponding optimally pure 
states. 

Example 1 Single qubit anisotropic depolarizing channel. 

Let T be the one-qubit channel given by p — > 
Ei=oP if7 ' P° l wnere the °" l ' s are the Pauli matrices (er° = 1) 
and the p^s a probability distribution. One finds S7qi — 
^io = y/PoPiV°v l = ^/PoPicr 1 ; Qij = -Qji = 
iy/PjPie tjk <J k (i = 1,2,3); flu = pjl, (i = 0, ...,3). It 
then follows that 

3 

fi(T)=^ai£7*<8>£7' (11) 
i=0 

where a = X)i=o^' a * = ^{PoPk + PiPj) (i + j + 
k,k = 1,2,3). Note that £f =0 Oj = (E?=o^) 2 = L that 
a G [1/4,1], and a fc e [0,l/2],(fc = 1,2,3). The eigen- 
states of Q(T) are the Bell states = (|01) - \10))/y/2 
(singlet) and {|<T) = (|00) - |11))/V5, |^+) = (|01) + 



\10))/V2, = (|00) + |H))/V2} (triplet). Their respec- 
tive eigenvalues are 2ao — 1 and {1 — 2a-y, 1 — 2«3, 1 — 20:2}. 
Further note that Specfi(T) C [-1/2, 1] and that fi(T) in 
the triplet sector is a positive operator. The triplet sector is 
symmetric, while the singlet is anti-symmetric. From Eq. @ 
we thus know that the minimal eigenvalue in the triplet sector 
provides a lower bound on the purity: 

P(T) > 1 - 2 max a % > (12) 

i=l,2,3 

(of course in this single-qubit example there are no non-trivial 
subspaces: C = fi). In this general case we cannot directly 
determine the actual purity or find the corresponding maxi- 
mally robust state(s), since the Bell triplet-states are not prod- 
uct states. To find the optimal purity states in such a case one 
has to resort to other optimization techniques. However, in 
certain special cases the eigenstates of il(T) will be product 
states, whence our method directly yields the optimally ro- 
bust states. For instance, consider the case po = p\ = 1/2 
and P2 — P3 = 0; one finds f2(T) = (1 + u x % cr x )/2. 
In this case \ip~),\<f)~) are degenerate, as are \ip + ), \<j> + ). 
We then find, respectively, the symmetric product eigenstates 
[(|0) - |l))/\/2p 2 and [(|0) + |l))/\/2]® 2 , both with eigen- 
value 1. The states |±) := (|0) ± \1))/V2 are thus both T- 
DFSs. This is intuitively clear, as the channel in this case is 
simply T(p) = p a p + p\cr x pa x ', and the states |±) lie on the 
Bloch sphere x-axis, which is invariant. 

As another example, consider the fully depolarizing chan- 
nel with pi = (1 — j?o)/3, (i — 1,2,3). Then the follow- 
ing (antiferromagnetic Heisenberg exchange) Hamiltonian is 
obtained: f2(T) = aol + a^ir' ® a 1 where ao = 
p 2 + (l- Po ) 2 /3,a = (2/3)[p (l-Po) + (1 ~Po) 2 /3] > 0. 
We can rewrite this as fi(T) = «ol + ct(2S — 1), where 
the SWAP operator S is defined by its action on basis states: 
S\i) (g> |i) = |i) 8> \i). In this case clearly every symmetric 
product state is an eigenstate of f2(T), with eigenvalue ao + a, 
which equals the channel purity. Thus all single-qubit states 
are equally (and optimally) robust. Again, this is intuitively 
clear: the fully depolarizing channel shrinks the Bloch sphere 
isotropically. 

Example 2 Correlated two-qubit anistropic depolarizing 
channel. 

Consider the correlated map 

T (p) = ^Pai^a <E> a a )p{a a (g) er Q ). (13) 
Then Q| a = ^/PaPpo'aO'p ® o- a crp and 

°= E pK 1 ® 1 )®^ E p a p f 3(<j^^r 2 . (H) 

a=0, x, y,z 

This example can be solved directly by observing that each of 
the Bell states has eigenvalue +1 or —1 under the action of 
<7 7 (8 <7 T , whence the purity is one. Thus the Bell states are 
T-DFSs. Notice that this result appears to be related to the 
communication problem for channels with correlated noise 
studied in Ref. Jivtl . 

One can also find the minimal purity states by differenti- 
ating (ft) :— ® 2 (V'|r2(T)|V>)' 82 as a function of expansion 
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parameters of \tp) over the Bell states. This yields (f2) m ; n = 
Sa=o x y z Pa> an< ^ me corresponding minimally robust set of 
states are superpositions of pairs of Bell states with arbitrary 
phases: 

|^) = — (|00) + |11)) + — (|01)-|10)) 

pia-z pifa 

l^ 2 > = — (|00)- |11)) + — (|01) + |10)). 

This channel thus has the interesting property that the maxi- 
mally entangled Bell states are more robust than any separable 
(pure) state. 

Example 3 Amplitude damping. 

Let T(p) = |0)(0|, Vp £ S(€ d ). A set of Ki-aus oper- 
ators is given by A, = |0)(i|, (i — l,...,d). Note that 
the channel is non-unital: X^AjA] = d|0)(0| > 1. One 
has A]Ai = \j){i\, so that Sl(T) = Eij ® KXj'I = 
= & ( the SWAP operator). Here all states 
are mapped onto a pure one and f2(T) is identically 1 in 
the symmetric subspace. A slight generalization is given by 
TO) = (l-p)p + p\0)(0\,(p£ [0,1])- In this case one finds 

Q(T) = (l-p)H+p 2 S 

+ p(l-j>)(|0)(0|<8l + l<8|0}<0|)£ (15) 

Note that ||fi(T)|| < (1 - p) 2 + p 2 + 2p(l - p) = 1. The only 
T-DFSisC|0). 

Example 4 Projective measurements. 

Let TO) = Ei n i/° n i. n * n j = *«n i ,X) i ni - 1- Then 
= AylL;, from which 

0(10=531^® H. (16) 

i 

If % := ImlL then Tif 2 C H Q , i.e., from Proposition 2 all 
the eigenvectors of the IT's are ID T-DFSs. The maximum 
eigenvalue of fi(T) is 1; this follows from (*|fi(T)|*) < 
1 (V|\E / )). The latter inequality results from the following ar- 
gument: 1® 2 = (E.IL)® 2 = IT ® IL, = fl(T) + 
Ylijtj n« <8> Ilj . The last term is a non-negative operator (sum 
of products of non-negative operators), whence 1 — fi(T) > 0. 
Taking the expectation value of the last inequality with respect 
to \^>) proves the bound above. 

In the following we use the operator norm ||A|| := 
max|||^)||| = i ||A|-i/>)||. We shall write \\A\\ for simplicity. 
Example 5 Unitary mixture of a group representation. 

This is a rather general and quite important example, which 
includes Examples 1,2 above. Let 

T{p) = Y J P 9 U gP Ul (17) 

9 

where g i— > p g is a probability distribution over the group 
C/ = {<?} and g i— > {/ g a unitary representation of . One finds 

= y/p^phWgUh = y/pJp~hU g -i h ; thus 

0(T) = ^g fc C/ fc OC/ fe -x, (18) 
keg 



where := J^gPgPgk- 1 is a l so a probability distribution. 
If |V>) e W is a ^-singlet, i.e., {7 9 |^>) = \ip) (Vc? G G) 
then |i/>)® 2 € ft n , i.e., all the ^-singlets are ID T-DFSs. 
Here again, the maximum eigenvalue of fi(T) is 1: Indeed, 
it is easy to see that || fi(T) 1 1 < 1: ||ft(T)|| < J2kVk\\U k ® 
U k -A\ = J2k1k = 1. 

As a particular instance of this kind of channel let us con- 
sider an A-qubit case with the Uk's generating an Abelian 
subgroup of Q the Pauli group (all tensor products of Pauli 
matrices on N qubits), as was the case in Examples 1,2 
above. The set of ^-singlets is now given by the the stabi- 
lizer of Q [denoted Stab(C/)], i.e., the subspace generated by 
the \ip) such that U k \ip) = \ip) (Vfe). Since U k = U\ one 
finds immediately that elements of the form \ip)® 2 , where 
|V-') € Stab((/), are eigenvectors of fi(T) with maximum 
eigenvalue ~^2 k q k . These states also play the role of code- 
words of stabilizer QECCs |8]. We thus see that, in this ex- 
ample, the stabilizer-QECC codewords are maximally robust, 
though no active error correction is assumed. 

In fact, the connection to quantum error correction can be 
made more general: the formalism developed so far allows 
us to establish an intriguing identity for the purity of states 
belonging to a QECC C for the CP-map T: p -> £\ A t pA\. 
If \tp a ): \ipp) S C then the error-correction condition is: 

(ipalAlAjltpp) = Cij8 a p, (19) 

where the matrix is Hermitian, non-negative, and has trace 
one 1 7]. For non-degenerate codes Cij has maximal rank. Let 
us now consider states of the form IV-'a)® 2 (IV't*) £ C). From 
Eq. and the error correction condition one has that 

p(t,c) = (^ 2 \n(T)\^ 2 ) = J2\(^\aUM\ 2 

ij 

= ^|cy| 2 = Tr(c 2 ). (20) 

Viewing c as a state (density operator), we have thus found 
that the purity of the channel acting on the codewords of C is 
just the purity of the "state " c associated to the code itself. 
For example, for DFSs c is simply a rank one matrix with unit 
trace 11411 . whence Trc 2 = 1 and the maximum eigenvalue 
condition is readily recovered. As a more interesting exam- 
ple, consider a CP-map T with A, = ^JpiUi with unitary C/j's 
(e.g., chosen from the Pauli group, as in stabilizer QEC). Re- 
call that c = A^A, where the matrix A is defined by the error- 
recovery relation R r Ai = X„t (restricted to C), for each re- 
covery operator R r 0|. Then R r — \riU~ 1 / \fpi, and from 

the CP-condition J2 r R\Rr = 1 we find J2 r \Ki\ 2 = Pi- As- 
sume for simplicity that there is a unique recovery operator 
per error, i.e., X r i — \i8 r i, Xi ^ Vi (this is an example of 
a non-degenerate code). Then |A;| 2 = pi and c = diagOi); 
it follows that the purity over such a QECC associated to T is 
simply given by ^p 2 . 



5 



V. THE DUAL REPRESENTATION 

We now develop an alternative representation of the 
channel-Hamiltonian, which is useful for the derivation of 
several additional results, and sheds new light on the physi- 
cal interpretation of the channel purity. 
Definition 6 The dual T* of a CP-map T [see Eq {!)] is: 
T*{X) = Y Ji A\XA i . 

Proposition 3 Let S be the SWAP operator (defined above). 
Then: 

VL{T)=Tf 2 (S)S. (21) 

We give two different proofs. 
Proof, a) 

T® 2 {S)S = Y^iA <8> A$)S(Ai <8> Aj)S 

ij 

= £(4 ®At)(A,-®A) 

ij 

= 53 A l A 3 ® A ] A * = n ( T ) w 

ij 

b) By writing the SWAP operator explicitly as S — 
J2 lm \m)(l\ ® |0( m l an d a PPly m g T® 2 one obtains: 

£im,y 4l m )('IA®^]lO(H^i- Then the proof follows by 
explicity comparing the matrix elements of the latter operator 
times 5 with the ones of ft(T). ■ 

We remark that one is led to consider the operator T® 2 (S 1 ) 
by the following argument: 

Tr[T 2 (p)] = Tj[S{T(p)®T(p)}] 

= Tr[ST® 2 (p®p)} = Ty{T® 2 (S)p® P }, 

where in the first step we used the identity 

Tr[AB] = Tr[5A ® B] (23) 

valid for general operators A, B fl5ll . and in the last step we 
"dualized" the map. Then for pure inputs p = \ip)(ip\ one has 
(V>® 2 |fi(T)|V>® 2 ) = (^® 2 \T® 2 (S*)NA® 2 )- This dualization 
is quite useful since it moves the burden of calculation of the 
channel action away from the entire set of states p to the single 
observable S. 

Coroiiary 1 Upon restriction to the to the symmetric subspace 
ofH® 2 one can write ft (T) = T® 2 {S). 
Proof. Immediate. ■ 

The following corollary contains a general derivation, 
based on the dual representation of Q(T), of a fact that was 
already proved for specific examples in Section lTVl 
Corollary 2 || ft(T) || < 1. 

Proof. One has || ft(T) || = ||T® 2 (S)S|| < ||Tf 2 (S)|| ||5|| < 
\\T® 2 {S)\\. Since T® 2 is the dual of a CP map elements 
smaller (greater) than the identity (minus the identity) are 
mapped onto elements smaller than the identity. Since — 1 < 
S < 1 one has -1 < T® 2 (S) < t. This relation implies 
in particular that the maximum eigenvalue of the Hermitian 
operator T® 2 (S*) is smaller than one. Since this maximum 



eigenvalue coincides with the IHI^ norm of T® 2 (S) the in- 
equality is proved. ■ 

We now present a result that allows one to directly compute 
the average purity of a quantum channel. 
Proposition 4 The Haar average purity of the CP-map T is 
given by 

W^WMf = d(dTT) TT[ST ® 2il) + n{T)] (24) 

Proof. Using the fact that J dip\ip) 2 is the normal- 
ized projector over the symmetric subspace of 7i® 2 , i.e., 
(t + S)/d(d+l) OJ], one has 

/ [Trmmn = 

= 1 Tr [ST® 2 (&)+T® 2 (S)S] (25) 



In other words, the Haar-average purity of a channel T is 
given by the expectation value of f2(T) over the maximally 
mixed state T\ + (7i) = l/d(d + 1) (1 + S) over the symmetric 
subspace ofH® 2 . 

Corollary 3 Using Eq. J24t one can get the Haar-averaged 
purities of the channels considered above: 

(i) One qubit depolarizing channel: (1 + 2ao)/3. 

(ii) Amplitude damping channel: (1— p) 2 +p 2 +2p(l—p) / d. 

(iii) Projective measurements: [d + X)i(Trn i ) 2 ]/d(rf + 1). 

(iv) Unitary mixing: [d + J2 ge g ig^U a?] / d ( d + l ) 
^From (iii) and (iv) it follows that: 

(a) One-dimensional projective measurements achieve the 
minimal average purity, of 2/{d + 1). 

(b) For unitary mixing and assuming a Haar-uniform distri- 
bution (all q g equal, i.e., the fully depolarizing channel), mini- 
mal purity is obtained for {7 ff 's in a tj-irrep. Indeed, one has in 
general that (1/|0|) Y, ge g |TrC/ | 2 = £j n 2 where nj is the 
multiplicity of the J-th Q irrep 1 16]. The minimum is clearly 
achived when just one irrep appears, i.e., the irreducible case. 

Before concluding this section we would like to point out 
that the formula (*|ft(T)|tf) = Tr [ST® 2 (|*)(*|)] allows 
us to give an operational meaning to the operator ft, and in 
the particular case in which |W) = \ip)® 2 , to the purity of 
T(\t/j)(ijj\). Indeed, this expectation value of ft(T) is noth- 
ing but the expectation value of the observable S in the state 
T® 2 (|*)(*|). The latter state can in turn be viewed as the 
result of an action of the channel on a pair of, possibly entan- 
gled, input states from Tt. 



VI. PURE STATE FIDELITY OF A CHANNEL 

We now show how many of the techniques introduced 
above for the channel purity carry over to the (simpler) prob- 
lem of calculating the 
Definition 7 Pure state fidelity: 

F(T,|^)) :=<V|I\WMM. (26) 
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Proposition 5 

(i) 



F(T, = 2 \n, (T) \^ 2 ) 



where fii(T) = (1 ® T*)(S)S. 
(ii) 



Jl 1 (T) = ^i 1 ®4. 



(27) 



(28) 



(iii) 



1 



d(d + l 



-Tr[Q 1 (T) + 5(l®T(l))]. (29) 



In particular, for a unital map one has average pure state fi- 
delity given by [d + £. |TrA ( | 2 ]/d(d + 1). 
Proof, 
(i) 

F(T,\i;)) = (^|r(|-0)<^|)|-0> = Tr [i^xvi^rdvX^D] 
= Tr [5|0) (0| (E) r(|^)(0|)] 

= Tr[(ll®T,)(5)|0)(0|» 2 ] 

= (^ 2 |(i®t;)(s)s|^ 2 ). 



(a) 



(1®T.)(5)5 = ®A\)S(l®Ai)S 

i 

= ^(i®4)(4®i) 



(iii) 



F(T,|^))" = / Tr[|V><V® 2 «i(T)] 

= Tr[l/d(d+l)(l + S)fii(T)] 
tr [ni(T) + (l®TQ(S)] 
d(d + l) 

Notice that the second term inside the square brackets is, 
for unital maps, simply TrS = d. 
■ 

Is important to stress that £!i (T) defined above is in general 
noM-Hermitian. On the other hand, Qi(T)S = (t ® T"*)(5) 
/s Hermitian (image of an Hermitian operator via CP-map) 
and has the same expectation values as f2i(T) over symmet- 
ric states in H® 2 . We thus associate to T a second channel 
Hamiltonian: 



Definition 8 The channel fidelity-Hamiltonian is: 

n'(T):=(l®r»)(5). (30) 

We now report, as corollaries of point (iii) above, the aver- 
age pure state fidelities of a few relevant channels. 
Corollary 4 

(i) Mixing of unitaries from the Pauli group (JV qubits): 
1/(2 N + 1)(1 + 2V 

(ii) Mixing of general unitaries: [d + ^ ■ pi\tiUi\ 2 } / d(d + 
1)). 

(iii) Amplitude damping: 1 — p(l — 1/d). 

(iv) Projective measurements: [d + Yli iTrlT^ | 2 ] / d(cZ + 1). 
As in the channel-purity case, the result i29\ can be simply 

stated by saying that: the Haar-average fidelity of a channel T 
is given by the expectation value ofQ\ (T) over the maximally 
mixed state IT + (Ti) = (1 + S)/ d(d + 1) over the symmetric 
subspace ofH® 2 . We note that a formula related to Eq. (I29> 
for the average fidelity of quantum operations has been given 
inRef. Gl. 



VII. CONCLUSIONS AND OUTLOOK 

We have introduced a Hamiltonian operator formalism for 
the calculation of the channel purity and pure state fidelity. 
Using this formalism we have been able to analytically com- 
pute these measures for a variety of channels of interest in the 
theory of open quantum systems, and quantum information 
theory. These analytical results are restricted to cases where 
the eigenstates of the Hamiltonian O (or ft') are product states 
in the symmetric subspace of H® 2 . When this is not the case 
one may have to resort to numerical methods to compute the 
purity/fidelity. 

A tempting generalization of our method is to consider 
perturbations to the channel Hamiltonian and use the well- 
developed tools of perturbation theory to thus study pertur- 
bations to given channels. One may further speculate about 
an adiabatic approximation, wherein slowly time-dependent 
channels can be studied using the adiabatic theorem applied 
to the channel Hamiltonian. We leave these as subjects for 
future investigations. 
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